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ABSTRACT 


Target recognition and identification in battlefields has been a crucial determinant 
to the ultimate success or failure of modern military campaigns. Since World War II, the 
Identification of Friend or Foe (IFF) systems installed in radar systems have served as the 
primary cooperative target identification techniques based on the “question and answer” 
interrogation loop of unidentified aircraft. However, the IFF system has a number of 
limitations that pose significant challenges in the positive identification of hostile and 
neutral aircrafts, which can lead to a catastrophic outcome of fratricide and the possible 
elimination of a friendly or commercial aircraft. 

To reduce the probability of fratricide and enhance the effectiveness and 
robustness of target identification, this research will examine the radar Non-Cooperative 
Target Recognition (NCTR) technique of using the bi-spectral signatures of backscattered 
radar signals. The basic idea is that the geometry of the target scatterers and their mutual 
interactions impose features in the reflected radar signal that are typical and unique to the 
target of interests. The bi-spectrum can be used to detect these multiple interactions 
features which may then be used to match against a reference database that contains 
signatures of different target types for recognition and identification. 


V 



THIS PAGE INTENTIONALLY LEET BLANK 


VI 



TABLE OF CONTENTS 


I. INTRODUCTION.I 

A. MOTIVATION.I 

1. IFF System and Cooperative Target Recognition.2 

2. Tragic Incidents of Fratricide.4 

3. Non-Cooperative Target Recognition (NCTR).5 

a. Electro-Optical (EO) and Infrared (IR) Systems . 5 

b. Electronic Warfare (EW) and Electronic Support 

Measure (ESM) Systems . 5 

c. Acoustics System . 6 

d. Radar System . 6 

B. SCOPE AND ORGANIZATION.7 

1. Thesis Scope.7 

2. Research Objectives.7 

3. Outline and Organization.7 

II. NCTR CLASSES AND TECHNIQUES.9 

A. FUNDAMENTAL RADAR TARGET RECOGNITION CLASSES.9 

B. HRR IMAGING.9 

1. Radial Range Profile.9 

2. 2D-ISAR Imaging.II 

C. SIGNAL MODULATION ANALYSIS.12 

I. JEM, PROM and HERM.12 

III. RADAR PRINCIPLES AND TARGET SCATTERING.15 

A. OVERVIEW OF RADAR CONCEPT.15 

1. Principle of Radar Operation.15 

2. Theory of Electromagnetic Wave Propagation.16 

3. Radar Signal.18 

B. TARGET SCATTERING MECHANISMS.19 

C. REPRESENTATION OF RADAR BACKSCATTERED SIGNAL.21 

1. Mathematical Model for Direct Specular Scattering.21 

2. Mathematical Model for Indirect Multiple Scatter Interactions ..22 

3. Mathematical Model for Combined Direct Specular and 

Indirect Multiple Scattering Interactions.23 

D. EXAMPLE: RADAR BACKSCATTERED SIGNAL FROM A 

SIMPLE 3-POINT SCATTERERS TARGET.24 

E. EXAMPLE: RADIAL RANGE PROFILE OF A SIMPLE 3-POINT 

SCATTERERS TARGET.25 

IV. BI-SPECTRAL NCTR TECHNIQUE.27 

A. OVERVIEW OF BI-SPECTAL METHOD.27 

1. Definition of Bi-Spectrum.27 

2. Definition of Bi-Time.28 

vii 






































3. Definition of Bi-Range.28 

B. EXAMPLE: BI-RANGE OF A SIMPLE 3-POINT SCATTERERS 

TARGET.29 

C. CHARACTERISTIC BI-RANGE PROFILE.32 

V. BI-SPECTRAL NCTR EXPERIMENTATION.35 

A. “BACKHOE DATA DOME” EXPERIMENTATION.35 

1. AFRL “Backhoe Data Dome Version I.O”.35 

2. Test Case I: Setup and Results (0azimuth=O°, (peievation=0°).37 

a. Radial Range Profile ( 0azimuth=O ° (peievation=0 °) . 38 

b. 3-D Bi-Range Profile (0azimuth=O ° (peievation=0 °) . 39 

c. 2-D Characteristic Bi-Range Profile (0admuth=O ° 

(Pelevation—0 ^. 41 

3. Test Case 2: Setup and Results (0azimuth=9O°, (peievation=0°).42 

a. Radial Range Profile ( 0azimuth=9O ° (Pelevation~ 0°) . 43 

b. 3-D Bi-Range Profile (0azimuth=9O ° <peievation=0 °) . 45 

c. 2-D Characteristic Bi-Range Profile (0azimuth=9O ° 

(Pelevation—0 °) . 46 

4. Test Case 3: Setup and Results (0azimuth=6O°, (peievation=30°).48 

a. Radial Range Profile ( 0azimuth=6O ° ^elevation =30°) . 49 

b. 3-D Bi-Range Profile (0azimuth=6O ° ^elevation~30 ^. 50 

c. 2-D Characteristic Bi-Range Profile (0azimuth=9O 

^elevation~ 0°) .52 

B. ANALYSIS OF “BACKHOE DATA DOME” EXPERIMENTATION ..53 

1. Radial Range Profiles.53 

2. Polarization Effects.55 

3. Characteristic Bi-Range Profile.55 

VI. CONCLUSIONS AND RECOMMENDATIONS.57 

A. SUMMING UP OF THESIS RESEARCH.57 

B. RECOMMENDATIONS FOR FUTURE WORK.58 

1. Experimentation Based on Actual Radar Signature 

Measurements of Real Complex Target.58 

2. Bi-spectral Waveform Design.59 

APPENDIX MATLAB SOURCE CODES.61 

A. SOURCE CODE FOR RADAR BACKSCATTERED SIGNAL.61 

B. SOURCE CODE FOR GENERATING THE THIRD-ORDER 

CUMULANT AND BI-SPECTRAL PLOTS.62 

C. SOURCE CODE FOR TEST CASE 1.63 

D. SOURCE CODE FOR TEST CASE 2.66 

E. SOURCE CODE FOR TEST CASE 3.68 

LIST OF REFERENCES.71 

INITIAL DISTRIBUTION LIST.73 

viii 





































LIST OF FIGURES 


Figure 1 IFF “Question and Answer” Loop (Wolff, 1997).2 

Figure 2 Radial Range Profile of a CAD modeled aireraft (Defenee Researeh and 

Development Canada).10 

Figure 3 ISAR Image of a B727 Jetliner (Borden, 1998).11 

Figure 4 JEM Speetrum of a Twin-Engine Aireraft (Van Der Heigen, 1998).13 

Eigure 5 PROM Speetrum of a DC-7 four-engine aireraft (Skolntk, 2001).13 

Eigure 6 HERM of a Hovering Helieopter (Bullard & Dowdy, 1991).14 

Eigure 7 Basie Radar Bloek Diagram and Operating Prineiple (Wolff, 1997).15 

Eigure 8 Electromagnetic Wave With Mutually Orthogonal E-Eield, H-Eield and 

Direction of Propagation Components.16 

Eigure 9 Target Scattering Mechanisms of Radar Wave (Tait, 2005).19 

Eigure 10 Huygens’ Principle of Diffraction (Potter).20 

Eigure 11 Propagation Path of Direct Specular Scattering Response from Individual 

Scatterers.21 

Eigure 12 Propagation Path of Indirect Scattering Response from Multiple 

Scatterers’ Interactions.22 

Eigure 13 Direct and Indirect Propagation Path of a 3-Point Scatterers Target.24 

Eigure 14 Radial Range Profile of 3-Point Scatterers Target.26 

Eigure 15 3-D Bi-Range Profile of 3-Point Scatterers Target.31 

Eigure 16 2-D Bi-Range Profile of 3-Point Scatterers Target.31 

Eigure 17 2-D Characteristic Bi-Range Profile of 3-Point Scatterers Target.32 

Eigure 18 CAD Model of Backhoe Eoader.35 

Eigure 19 Backhoe Data Dome Representation in k-Space (Ein & Naidu, 2004).36 

Eigure 20 Orientation of CAD-Modeled Backhoe Loader at Aspect Viewing Angles 

of Azimuth, (0azimuth)=O° and Elevation, ((peievation)=0°.37 

Eigure 21 Radial Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=O° 

and (peievation=0° with VV Polarization.38 

Eigure 22 Radial Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=O° 

and (Peievation=0° with HH Polarization.38 

Eigure 23 Radial Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=O° 

and (peievation=0° with HV Polarization.39 

Eigure 24 3-D Bi-Range Profde of Backhoe Eoader at Aspect Angles of 0azimuth=O° 

and (peievation=0° with VV Polarization.39 

Eigure 25 3-D Bi-Range Profde of Backhoe Eoader at Aspect Angles of 0azimuth=O° 

and (peievation=0° with HH Polarization.40 

Eigure 26 3-D Bi-Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=O° 

and (peievation=0° with HV Polarization.40 

Eigure 27 2-D Characteristic Bi-Range Profile of Backhoe Eoader at Aspect Angles 

of 0azimuth=O° and (peievation=0° with VV Polarization.41 

Eigure 28 2-D Characteristic Bi-Range Profile of Backhoe Eoader at Aspect Angles 

of 0azimuth=O° and (peievation=0° with HH Polarization.41 


IX 































Figure 29 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect Angles 

of 0azimuth=O° and (peievation=0° with HV Polarization.42 

Figure 30 Orientation of CAD-Modeled Backhoe Loader at Aspect Viewing Angles 

of Azimuth, (0azimuth)=9O° and Elevation, ((peievation)=0°.42 

Figure 31 Radial Range Profile of Backhoe Loader at Aspect Angles of 0azimuth=9O° 

and (Peievation=0° with VV Polarization.43 

Figure 32 Radial Range Profile of Backhoe Loader at Aspect Angles of 0azimuth=9O° 

and (peievation=0° with HH Polarization.44 

Figure 33 Radial Range Profile of Backhoe Loader at Aspect Angles of 0azimuth=9O° 

and (Peievation=0° with HV Polarization.44 

Figure 34 3-D Bi-Range Profde of Backhoe Loader at Aspect Angles of 0azimuth=9O° 

and (peievation=0° with W Polarization.45 

Figure 35 3-D Bi-Range Profde of Backhoe Loader at Aspect Angles of 0azimuth=9O° 

and (Peievation=0° with HH Polarization.45 

Figure 36 3-D Bi-Range Profile of Backhoe Loader at Aspect Angles of 0azimuth=9O° 

and (peievation=0° with HV Polarization.46 

Figure 37 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect Angles 

of 0azimuth=9O° and (peievation=0° with VV Polafization.46 

Figure 38 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect Angles 

of 0azimuth=9O° and (peievation=0° with HH Polafization.47 

Figure 39 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect Angles 

of 0azimuth=9O° and (peievation=0° with HV Polafization.47 

Figure 40 Orientation of CAD-Modeled Backhoe Loader at Aspect Viewing Angles 

of Azimuth, (0azimuth)=6O° and Elevation, ((peievation)=30°.48 

Eigure 41 Radial Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=6O° 

and (Peievation=30° with W Polarization.49 

Eigure 42 Radial Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=6O° 

and (peievation=30° with HH Polarization.49 

Eigure 43 Radial Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=6O° 

and (Peievation=30° with HV Polarization.50 

Eigure 44 3-D Bi-Range Profile of Backhoe Loader at Aspect Angles of 0azimuth=6O° 

and (peievation=30° with VV Polarization.50 

Eigure 45 3-D Bi-Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=6O° 

and (Peievation=30° with HH Polarization.51 

Eigure 46 3-D Bi-Range Profile of Backhoe Eoader at Aspect Angles of 0azimuth=6O° 

and (peievation=30° with HV Polarization.51 

Eigure 47 2-D Characteristic Bi-Range Profile of Backhoe Eoader at Aspect Angles 

of 0azimuth=6O° and (Peievation=30° with VV Polafization.52 

Eigure 48 2-D Characteristic Bi-Range Profile of Backhoe Eoader at Aspect Angles 

of 0azimuth=6O° and (peievation=30° with HH Polafization.52 

Eigure 49 2-D Characteristic Bi-Range Profile of Backhoe Eoader at Aspect Angles 

of 0azimuth=6O° and (peievation=30° with HV Polafization.53 

Eigure 50 Comparison of Radar Range Profile with Actual Backhoe Eoader.54 


X 

























LIST OF TABLES 


Table 1 Summary of Radial Range Profile’s Advantages and Disadvantages.10 

Table 2 Summary of 2D-ISAR Imaging’s Advantages and Disadvantages.12 

Table 3 Summary of JEM, PROM and HERM Analysis’s Advantages and 

Disadvantages.14 

Table 4 Summary of Bi-range Profile’s Advantages and Disadvantages.33 







THIS PAGE INTENTIONALLY LEET BLANK 



ACKNOWLEDGMENTS 


I wish to thank my wife, Charlene, for her love and eneouragement throughout the 
eourse of my researeh and study at the NPS. 

I would like to thank my thesis advisor. Professor Brett Borden, for his exeellent 
guidanee and knowledge-sharing during my entire thesis research and learning process. I 
also wish to thank my thesis co-advisor. Professor Donald L. Walters, for his invaluable 
advice and discussions. 

I am thankful to my fellow coursemates. Major Charles Chua, Major Anthony 
Lee, Captain James Neo and Mr. Yeo Eng Choon, for their friendship, encouragement 
and support. 



THIS PAGE INTENTIONALLY LEET BLANK 


XIV 



I. 


INTRODUCTION 


Know the enemy and know yourself; in a hundred battles you will never 
be in peril. Know yourself but not the enemy; for every vietory gained 
you will also suffer a defeat. Know neither the enemy nor yourself; you 
will sueeumb in every battle. 

— Sun Tzu, Chinese General and Military Strategist (1971) 


A. MOTIVATION 

Situation awareness has always been the key determinant for the ultimate victory 
or defeat of all warfare and military campaigns. Through situation awareness, one must 
be able to ascertain the four “Ws” of one’s own forces as well as one’s enemies, which 
are: 

1. Where are they: azimuth, elevation, range, etc? 

2. When will they arrive: heading, velocity, time of flight, etc? 

3. Who are they: friendly, hostile or neutral? and 

4. What are they: fighter jet, commercial aircraft, tank, school bus, etc? 

While “Where” and “When” are adequately addressed by the early RAdio 
Detection And Ranging (Radar) systems used primarily as sensors to detect and track 
targets since World War II, “Who” and “What” or more specifically, target recognition 
and identification, can pose significant challenges, as only limited information pertaining 
to the characteristics of the target being detected has been available from the radar 
backscattered signal. 

As such, early radar systems were merely “blob” detectors, as detections of target 
presences from the reflected radar signal were usually displayed as phosphor “blips” on 
the radar system’s Plan Position Indicator (PPI) (Skolnik, 2001). Information on the 
target’s location in range, azimuth and elevation and velocity could then be extracted 
from the PPI, but could not provide further detail on the target’s type and identity. 
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This limitation of the early radar systems, together with the vital need for situation 
awareness as a foree multiplier to enhance combat effectiveness, provided the impetus for 
the development of the IFF system and various other radar target recognition techniques. 

1. IFF System and Cooperative Target Recognition 

The IFF system is a two-channel interrogator / transponder system, with one 
frequency used for the interrogating signals and another for the response, that uses pulse- 
coded radio frequency transmissions to positively identify friendly aircraft, ships, and 
army units from those of the enemy forces. It serves as a cooperative technique for target 
recognition, based on the “Question” interrogation and “ Ans wer” response loop of 
unidentified targets of interests, as depicted in Figure 1. 



Figure 1 IFF “Question and Answer” Loop (Wolff, 1997) 


The IFF interrogator system initiates the “Question and Answer” loop by 
transmitting a pulse-coded radio frequency signal, also referred to as the “challenge” that 
asks a series of questions, one of which is “Who are you?” to the unidentified target of 
interest. The interrogation code or challenge is then received by an electronic system 
known as a transponder that is installed onboard the target. If the target is friendly and 
equipped with a compatible transponder, it automatically interprets the challenge and 
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responds by transmitting the requested identification code back to the interrogator 
system. After successfully receiving and decoding the acknowledgement, the target is 
tagged as a friendly target by the interrogating radar system. 

Conversely, if the target is neutral or hostile, it is assumed that the challenge 
cannot be decoded and interpreted by the onboard transponder. Henceforth, the 
transponder is unable to respond with a correct answer and will be tagged as an unknown 
or foe target by the interrogating radar system. 

Although the IFF system plays a significant role in contributing toward target 
recognition, it can only perform positive recognition of friendly targets, and it requires 
the active and cooperative participation of the targets which are under interrogation. 

Moreover, IFF systems have a number of limitations, and cooperative systems 
have historically been plagued with the following fundamental issues: 

a. If the IFF system is inoperative, due to a malfunction, or an operator 
intentionally turns it off, there is no inherent method for a weapon system to 
distinguish the resulting lack of response from that of an enemy target. 

b. Hostile targets may intercept / listen to and re-transmit an 
acknowledgement from a “friendly” target, in which case the hostile target will be 
tagged as “friendly.” 

c. Commercial airliners are fitted with transponders that conform to open 
international civil aviation standards, which enable terrorists or enemy forces the 
opportunity to be equipped with such transponders on their aircrafts. 
Consequently, they can masquerade as passenger airliners for the purpose of 
infiltrating the air defence system. 

Hence, it is obvious that the IFF system is not infallible and any interruption in 
the IFF “Question and Answer” loop, resulting in the absence of a response or correct 
acknowledgement, can potentially draw the inference that the unidentified target is 
hostile. This seriously reduces the utility of neutral or allied forces and can potentially 
lead to catastrophic outcomes of fratricide and the erroneous elimination of a friendly or 
commercial aircraft. 


3 



2. Tragic Incidents of Fratricide 

To illustrate the devastating consequences of erroneous target recognition and 
identification, a few examples of wrongful engagements of neutral civilian and friendly 
targets from recent history are referenced from V.G Nebabin’s Methods and Techniques 
of Radar Recognition (1995) and P. Tail’s Introduction to Radar Target Recognition 
(2005). 

a. In 1980, a DC-9 Italian passenger plane was shot down by mistake with a 
missile (supposedly American or French made) near the Island of Sicily. It is 
believed that the missile was re-aimed by mistake because of incorrect IFF and 
discrimination of the type of plane. Eighty-one passengers died. 

b. In 1983, a Boeing 747 passenger plane. South Korean Airlines Flight 007, 
flying from New York and Anchorage to Korea, was shot down by the USSR Air 
Defense SU-15 fighter as a result of incorrect recognition. The Boeing 747 was 
recognized as an RC-135 reconnaissance aircraft. Two-hundred and sixty-nine 
passengers from Korea, Japan, the U.S. and other countries died. 

c. In 1988, an Iranian A-300B Airbus passenger plane. Flight 655, flying 
from Bender-Abba, Iran, to Dubai, United Arab Emirates was shot down by an 
SM-2 Standard Missile under control of the Aegis system of the cruiser USS 
Vincennes in the Persian Gulf The A-300B was recognized erroneously as an E- 
14 fighter. Two-hundred and ninety-eight people died. 

d. In 1994, two U.S. Blackhawk military helicopters were shot down in the 
north of Iraq by two U.S. E-I5 fighters. The passengers and crews of the 
helicopters had been performing a humanitarian mission. They had been 
recognized erroneously as Iraqi helicopters flying in the no-fly zone. Eourteen 
men died. 

e. During the 2003 Iraq war, there were several incidents of ground troops 
being killed by friendly fire from the air. A UK Tornado aircraft was shot down 
by a Patriot missile, apparently being mistaken for a Scud ballistic missile, and 
two British soldiers were killed apparently in an exchange of fire between two 
British Challenger tanks. 
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3. Non-Cooperative Target Recognition (NCTR) 

To prevent the eatastrophie outeome of fratrieide and enhanee the robustness of 
the air defense system, NCTR, whieh eneompasses the teehniques and teehnologies to 
extraet and examine the eharaeteristies of target signatures in order to provide further 
insight on the target’s type and identity without its aetive eooperative partieipation, has 
been developed and adopted over the years (Tait, 2005). 

A variety of aetive and passive sensor systems ean be utilized to perform NCTR 
and establish a target’s type and identity. A diseussion of the various sensor systems and 
their strengths and limitations in target reeognition performanee is presented in the 
following seetions. 

a. Electro-Optical (EO) and Infrared (IR) Systems 

An EO system uses opties for the surveillanee, traeking and reeognition of 
a target and depends on external illumination for the eonversion of photons into amplifier 
eleetrieal signals for imaging. An IR system deteets the heat emission of the target and 
ean operate in total darkness without external illumination of the target as it is dependent 
on the temperature differenees between the target and its baekground, for imaging. Both 
EO and IR systems are passive sensor systems that are ideal for eovert operations and 
have strength in providing high-resolution images for target reeognition. However, their 
performanees are adversely affeeted by weather effeets, sueh as rain, fog and haze, and 
by atmospherie attenuation due to seattering and absorption. 

b. Electronic Warfare (EW) and Electronic Support Measure 
(ESM) Systems 

EW and ESM systems use passive reeeivers for the deteetion and analysis 
of signal in the eleetromagnetie speetrum. The eharaeteristies of the intereept signal are 
measured and used to isolate, identify and loeate the emitter. Identifieation of the emitter 
through the use of finger-printing teehniques of its measured signal, along with a library 
of referenee signals ean then be used to reeognize the assoeiate platform. Although EW 
and ESM systems have good range performanee, due to the one-way propagation of the 
emitter signal, their effeetiveness is eompletely reliant on the target’s eleetromagnetie 
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emission for locating and possibly identifying it. Hence an Emission Control (EMCON) 
imposed by an incoming target will easily overcome the EW and ESM systems’ NCTR 
capability. 

c. Acoustics System 

The acoustics system uses an array of microphones to discern sound 
waves from the noisy environment for the detection, localization and possible recognition 
of sound generating targets, such as human speech, the engine of a moving vehicle and 
the chopping blades of a hovering helicopter, etc. While the acoustics signal from these 
targets may be discriminating and recognisable, the acoustics sensor is generally limited 
in performance due to the slow speed of sound, which is about 330 m/s in air, and limited 
sensitivity to background noise. The acoustics system is also affected by atmospheric 
propagation attenuation and is only considered for very short-range applications. 

d. Radar System 

By far, a radar system is considered the sensor of choice for NCTR, as it 
offers the most potential. It is able to perform round-the-clock 24/7 surveillance, can 
perform detection and tracking of targets at ranges much further than EO and IR systems, 
and is not affected by adverse weather conditions. Eurthermore, unlike EW and ESW 
systems, radar systems are not dependent on the target emissions. Technological 
advancements have created a system that formerly could only detect and track targets to 
one that is capable of high resolution imaging and target recognition. 

According to Van Der Heigen’s “Aircraft Recognition with Radar Range 
Profiles (1998),” Radar NCTR techniques are rooted on the basic concept that the 
geometry of the target or the moving parts on the target impose characteristic features in 
the reflected radar signal that are typical and are unique to the targets of interest. These 
features may then be used to match against a reference database that contains signatures 
of different target types for target recognition and identification. 
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B. SCOPE AND ORGANIZATION 

1. Thesis Scope 

This thesis first discusses the two fundamental methods of High Range Resolution 
(HRR) imaging and Signal Modulation analysis, developed to perform target recognition 
and identification using a radar system. Several conventional NCTR techniques based on 
these two fundamental methods are then reviewed and evaluated based on their relative 
merits and limitations. 

Next, the thesis focuses on developing a complementary radar NCTR technique of 
using the bi-spectral signatures of backscattered radar signals to augment the knowledge 
that conventional NCTR techniques already offer. Derivation and theoretical calculations 
of the bi-spectrum of radar target signature were performed to extract the multiple 
scattering interaction features which may then be used to match against a reference 
database that contains signatures of different target types for target recognition and 
identification. 

2. Research Objectives 

The objectives of this thesis are two-fold: 

a. To study and mathematically model the direct and indirect multiple 
interactions scattering mechanisms of radar backscattered signal. 

b. To develop the radar NCTR technique of using the bi-spectral 
signatures of backscattered radar signals for characterization of a target’s multiple 
interactions features. 

3. Outline and Organization 

Chapter II provides a review of the current conventional NCTR techniques based 
on the two fundamental classes of HRR imaging and Signal Modulation analysis, 
developed to perform target recognition and identification using radar systems. A 
comparison of the NCTR techniques and their relative advantages and disadvantages in 
target recognition performance is presented. 
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Chapter III examines the radar principles and scattering mechanisms of 
electromagnetic waves when a radar signal intersects any target in its propagation path. 
A mathematical model to represent direct specular and indirect multiple scattering 
interactions of the backscattered radar signal is presented. 

Chapter IV introduces the definitions and mathematical derivations of the bi¬ 
spectrum, bi-time and bi-range of the radar backscattered signal. The concept of using the 
bi-spectral signatures, also known as the characteristic bi-range profiles of backscattered 
radar signal for target recognition and identification, is presented and studied through an 
example of a simple 3-points scatterer target. 

Chapter V further investigates the potential of using bi-spectra technique for radar 
target recognition by conducting experimentation based on the Air Force Research 
Laboratory, “Backhoe Data Dome.” The results and findings for three different test cases 
are presented and analyzed. 

Chapter VI discusses the conclusions made in this thesis and recommends 
objectives for follow-on research work. 
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II. NCTR CLASSES AND TECHNIQUES 


A. FUNDAMENTAL RADAR TARGET RECOGNITION CLASSES 

NCTR is generally eategorized into two fundamental elasses: High Range 
Resolution (HRR) imaging and Signal Modulation analysis to perform radar target 
reeognition and identifieation (Tait, 2005; Van Der Heigen, 1998). The HRR imaging 
elass of NCTR exploits the radar signatures of a target’s physieal shape, while the Signal 
Modulation analysis elass of NCTR exploits the radar signatures of the target or the 
dynamies of its moving-parts. A diseussion of the eoneepts, relative merits and 
limitations of the eurrent eonventional NCTR teehniques based on these two fundamental 
elasses is presented in the following seetions. 

B. HRR IMAGING 

The basis of the HRR imaging elass of NCTR is established on the ability for high 
radial and eross range resolutions radar to resolve the individual seattering eenters of a 
target and perform radar imaging for target reeognition. NCTR teehniques from this 
family inelude one-dimensional HRR profile or, simply, radial range profile and two- 
dimensional Inverse Synthetie Aperture Radar (2D-ISAR) imaging (Skolnik, 2001; Tait, 
2005; Van Der Heigen, 1998). 

I. Radial Range Profile 

The radial range profiles of a target of interest essentially represent the one¬ 
dimensional images of the target as viewed from different aspeet angles. The dominant 
individual scatterering eenters or subeomponents, i.e. the parts on the target that present 
large Radar Cross-Seetions (RCS) and subsequently return a strong radar refieetion, are 
projeeted onto the line-of-sight between the radar system and the target. 

An example of a one-dimensional image that ean be ean be ereated by the HRR 
radar system is illustrated in Figure 2. Figure 2 shows the radial range profile of a target 
generated by the eomputer predietion eode using the Computer Aided Design (CAD) 
model of an aireraft viewed from the left hand side. Impulse responses, ereated using the 
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Fourier transform of the backscattered radar signal from the aircraft dominant scatterers 
(dots), are then projected onto the line-of-sight in order to generate a radial range profile. 



Figure 2 Radial Range Profile of a CAD modeled aircraft (Defence Research 

and Development Canada). 


The radial range profile contains information on the geometry and heuristic 
features of the target and can be used to infer the type characteristic of the target. This 
information, together with collected reference databases of measured and / or stimulated 
(i.e. generated using RCS prediction tools and CAD models) target signatures of different 
target types (e.g. aircraft, ship and vehicle) and orientations, can then be processed by 
target recognition and identification algorithms to classify the unknown target into 
broader categories, such as small jet, large jumbo jet, missile and tank, etc. down to the 
specific platform type, such as F-16 fighter jet, Boeing-747 passenger plane, Kh-31 Anti- 
Radiation Missile and M-1 Main Battle Tank, etc. 

The merits and limitations of the radial range profile NCTR technique according 
to Van der Heigen (1998)are summarized in Table 1; 


Advantages 

Disadvantages 

a. Applicable to all target aspect angles. 

a. Radial range profile is highly dependent on 

the target aspect angle. Hence large data sets of 

target signatures covering all expected target types 

and orientations are required. 

b. Requires relatively short time on target. 


Table 1 Summary of Radial Range Profile’s Advantages and Disadvantages 
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2. 2D-ISAR Imaging 

Unlike the radial range profile which only uses the line-of-sight or down range 
information of the backscattered radar signal for HRR imaging, 2D-ISAR imaging uses 
information in both the down range and cross-range directions to form two-dimensional 
images of the target of interests. The cross-range resolution in 2D-ISAR can be thought 
of in terms of the high resolution in the Doppler-frequency domain. Since each of the 
major individual scatterer of a moving target can have a different relative velocity with 
respect to the radar system, resolution in Doppler-frequency, especially for the rotation 
components of the target’s motion, will make it possible to separate these scatterers in the 
cross-range direction. Figure 3 shows an example of an ISAR image. The target is a B727 
jetliner with orientation as in inset (Borden, 1998). 



crofss-rEmgCj meters 

Figure 3 ISAR Image of a B727 Jetliner (Borden, 1998). 
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The merits and limitations of the 2D-ISAR imaging NCTR teehnique according to 
Van Der Heigen (1998) are summarized in Table 2: 


Advantages 

Disadvantages 

a. Better fidelity in mapping the geometry of 

the target due to the additional resolution in the 

cross-range dimension. 

a. Requires rotational components of target’s 

motion perpendicular to the line-of-sight. Hence 

target with weak angular rotation (i.e. pitch, roll and 

yaw motions) will have poor cross-range resolution. 

b. More photographic-like 2D image which is 

better suited for human interpretation. 

b. Requires complex motion compensation 

and auto-focusing algorithm of 2D-ISAR images 

which is computationally heavy. 

c. Limited knowledge on the precise project 

plane for the rotation vectors, making recognition 

complicated. 


Table 2 Summary of 2D-ISAR Imaging’s Advantages and Disadvantages 


C. SIGNAL MODULATION ANALYSIS 

The Signal Modulation analysis class of NCTR is based on the characterizations 
of the reflected radar modulation signals from the rotating parts of the engine, such as the 
compressor blades in a jet engine, propellers on a propeller-aircraft and rotors on a 
helicopter, in the Doppler domain. The NCTR techniques from this family are 
prominently used to perform target recognition and the identification of aircrafts and 
include the following: Jet Engine Modulation (JEM), Propeller ROtor Modulation 
(PROM) and HElicopter Rotor Modulation (HERM) (Skolnik, 2001; Tait, 2005; Van Der 
Heigen, 1998). 


I. JEM, PROM and HERM 

The compressor blades in a jet engine, aircraft propellers and rotating blades and 
hub of a helicopter, each produces a unique and distinctive modulation of the radar target 
signatures. These modulations in the characteristic spectrums of JEM, PROM and 
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HERM can be analysed, and attributes pertaining to the information on engines / 
propellers / rotors types, quantities, number of blades, ranges of rotation frequencies, ete, 
can be extracted. Examples of the JEM speetrum of a twin-engine aircraft, PROM 
spectrum of a DC-7 four-engine aireraft and HERM speetrum of a hovering helieopter 
are shown in Eigures 4, 5 and 6 respeetively. 



Eigure 4 JEM Speetrum of a Twin-Engine Aireraft (Van Der Heigen, 1998) 



Doppler frequency (Hz) 

Eigure 5 PROM Spectrum of a DC-7 four-engine aircraft (Skolnik, 2001) 
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Figure 6 HERM of a Hovering Helicopter (Bullard & Dowdy, 1991) 


Since each aircraft has a unique engine, propeller or rotor, these characteristic 
parameters can be matched against a reference database or table that contains 
specifications for a variety of target engines / propellers / rotor types for engines 
associated target recognition and identification. 

The merits and limitations of JEM, PROM and HERM analysis NCTR techniques 
according to Van Der Heigen (1998) are summarized in Table 3: 


Advantages 

Disadvantages 

a. Small and simple target database since it is 

comprised of engine specification rather than radar 

target signature or image. 

a. Dependent on the target aspect angle. 

Useful spectrums mostly in the head-on or tail-on 

orientation. 

b. Requires relatively short time on target. 

b. Requires strong signal-to-noise ratio for 

effective classification. Hence techniques limited to 

relatively short ranges 

c. Fast classification since almost all aircrafts 

have a unique engine type. 


Table 3 Summary of JEM, PROM and HERM Analysis’s Advantages and Disadvantages 
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III. RADAR PRINCIPLES AND TARGET SCATTERING 


A. OVERVIEW OF RADAR CONCEPT 

RAdio Detection And Ranging (RADAR) is a method of using electromagnetic 
waves to remotely sense the presence and extract information pertaining to the position, 
velocity and identifying characteristics of targets of interest. An overview of the basic 
radar operating principle and electromagnetic wave propagation theory is presented in the 
following sections. 

1. Principle of Radar Operation 

The basic radar block diagram and principle of radar operation is illustrated in 



transmitter-^ duplexer-^ antenna electromagnetic wave-^aim 

display receiver^ duplexer ♦ antenna echo signal aim 

Figure 7 Basic Radar Block Diagram and Operating Principle (Wolff, 1997) 

Using a transmitter, the radar system first generates an electromagnetic signal 
(such as a short sinusoidal pulse) that is radiated into free space through a duplexer and 
an antenna. The duplexer switches the antenna of monostatic single-antenna systems 
between the transmitter and receiver, allowing the antenna to be shared for both 
transmitting and receiving functions, while the antenna is used to concentrate and direct 
the electromagnetic energy towards the target. 
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As the electromagnetic wave propagates in free space, a fraction of the 
transmitted electromagnetic energy is intercepted by the target and scattered in all 
directions. The reflected energy of the echo signal directed back towards the radar 
system, which is also known as the radar backscattered signal, is then collected by the 
same antenna and redirected to the receiver through the duplexer. Thereafter, the 
receiver performs the signal processing of the received echo signal in order to detect the 
presence of the target and extract information pertaining to its position, velocity and 
identifying characteristics. The process is completed with the presentation of this 
information on an operator display. 

2. Theory of Electromagnetic Wave Propagation 

The fundamental laws of electricity and magnetism, also known as Maxwell’s 
equations, predicted the phenomenon of electromagnetic waves in the form of sinusoidal 
electric E-field and magnetic H-field, which are oscillating at the wave’s frequency and 
mutually orthogonal to the direction of propagation, as depicted in Figure 8. 

Direction of 



Figure 8 Electromagnetic Wave With Mutually Orthogonal E-Field, H-Field 

and Direction of Propagation Components 
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By convention, the electrie E-field is usually adopted in the representation of the 
eleetromagnetic wave. A derivation of the eleetric E-field wave equation from 
Maxwell’s equations is presented according to Sadiku’s Elements of Electromagnetics 
(2001). Beginning from Maxwell’s equations for a vaeuum: 


v-e; = 0 

dt 

VH =0 

d 

VxH=£.—E 

°dt 


where 

V is the vector differential operator, 

jUq = 47t: X 10'^ H/m is the permeability of the vacuum and 
12 

£o = 8.854 X 10' E/m is the permittivity of the vacuum. 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


Taking the eurl of equation (3.2): 


Vx(VxE:) = V; 


d rj 

""“a" 


(3.5) 


Evaluating the left hand side of equation (3.5) with the vector identity of 
V X (v X a) = V(V • a)-V^A and simplifying the expression using equation (3.1), 


vx(vxe:) = v(v-e:)-v^e: = -v^£: 


(3.6) 


Evaluating the right hand side of equation (3.5) and simplifying the expression using 
equation (3.4), 


Vx 


d rr 


a 


a' 




(3.7) 
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Equating equations (3.6) and (3.7), results in a differential equation for the electric E- 
field wave equation, 

2 52 

S7^E-^oSo—E=0 (3.8) 

dt^ 

where 

jUgSo =\lc and 

c = 2.998 X 10^ m/sec is the speed of light. 

3. Radar Signal 

By solving the electric E-field wave equation in equation (3.8), the radar signal 
for an electromagnetic wave propagating from a radar system can be represented as a 
general form of a complex sinusoid electric E-field varying in both time and space and 
expressed as 

E{r,t) = A{co\cos{cot -k ■r)-i sin{cot - A: • r)] (3.9) 

where 

A is the amplitude of the E-field signal, 

CO =27rfis the radial frequency, 

k = — is the wave number, 

X is the wavelength, 

r is the range from the radar system and 
t is the time duration elapsed. 

Using Euler’s formula, = cosx - i sinx , the radar signal at a given position and time, 
E{r, t) can be rewritten in the exponential form as 

E{r,t) = A{(o)e~^'^^^'^'''^ (3.10) 
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To compute the eeho signal collected at the radar receiver, the range from the radar 
system is set to zero, r= 0, hence the received radar signal in the time domain is 

= (3.11) 

While the corresponding reeeived radar signal in the frequency domain is 

(3.12) 

B. TARGET SCATTERING MECHANISMS 

The mechanisms by which target scattering occurs, are important considerations 
for determining the physical size, shape, orientation and material composition of the 
target, which then account for the portion of electromagnetic energy that is reflected back 
in the direction of the radar system and consequently the probability of target detection, 
recognition and identification. 

As discussed in Tait (2005), the main interactions of the electromagnetic wave 
with the target are; reflection; absorption; refraction; diffraction and creeping surface 
waves. Figure 9 illustrates the interaction effects of reflection, absorption and refraction. 
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At the surface boundary of the target, the incident electromagnetic energy is 
reflected, refracted and absorbed by the interactions of the radar wave with the target. It 
can be observed that the fraction of incident electromagnetic energy reflected back 
towards the radar system is dependent on the portion of energy that is absorbed and 
refracted into the material, as well as the angle of incidence to the material surface. 

Another form of scattering mechanism, through the interaction of the radar wave 
with the target, is diffraction. This knife-edge effect is explained by Huygens' Principle, 
which states that a well-defined obstruction to an electromagnetic wave acts as a 
secondary source and creates a new wavefront. This new wavefront then propagates into 
the geometric optics shadow area of the obstacle, enabling the radar wave to bend around 
the edges of the target, propagate into the target’s structure through openings and detect 
the interior parts of the target that would be considered as non line-of-sight as shown in 
Figure 10. 



Figure 10 Huygens’ Principle of Diffraction (Potter) 
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Finally the creeping surface wave interaction effect allows for the propagation of 
radar electromagnetic wave along the surface of a target using surface wave to access the 
geometric optics’ blind region. 

C. REPRESENTATION OF RADAR BACKSCATTERED SIGNAL 

Targets are generally composed of many individual discrete specular scatterers, 
each with different scattering properties which correspond to a physical feature of the 
target, for example the nose, engines, main and tail wings of an aircraft. 

Each scatterer, when intercepted by a radar electromagnetic wave, will reflect an 
echo signal characterized by an amplitude and a phase back towards the direction of the 
radar system, which will then sum up to a resultant radar backscattered signal at the radar 
receiver. In addition, there will be contributions from the interactions between the 
multiple individual scatterers that will further complicate the target characteristics. 

A simple mathematical model for the direct specular and indirect multiple 
scattering interactions of the backscattered radar signal is discussed and presented in the 
following sections as referenced from Borden (1998) and Jouny (1994). 

I. Mathematical Model for Direct Specular Scattering 

The radar backscattered signal is assumed to be composed of two major 
components. The first component represents the specular direct scattering response from 
the individual scatterers, distributed over the target. The propagation path of the signal 
can be represented in Figure 11 as; 



tp/2 one-way propagation time 



pth seatter point 


tp/2 one-way propagation time 


Radar 


Round-trip propagation time = tp/2 -I- tp/2 = tp 


Figure 11 Propagation Path of Direct Specular Scattering Response from 

Individual Scatterers 
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Using equation 3.12, this component of the backscatter signal can be expressed 
mathematically as a summation of the reflected echo signals from all of the individual 
discrete specular scatterers; 


where 


hM= 1 . 4 , 

p=\ 


(3.13) 


m is the number of scattering centers, 

Ap is the amplitude of the scattered signal from the /?th scattering centre, 

y/p is the phase shift due to /?th scatterer and statistically characterized by the 
uniform probability density distribution over the interval (0, 27r), 

CO is the radial frequency {=27f) and 

tp is the round-trip propagation time between the radar and the />th scattering 
centre. 


2. Mathematical Model for Indirect Multiple Scatter Interactions 

The second major component of the radar backscattered signal represents the 
effect of interaction between multiple scattering points. The U* order interaction between 
any two scattering points is represented in Figure 12 as: 


ps 

Radar 


tp/2 one-way propagation time 


tpq interaction time 


pth scatter point 


tq/2 one-way propagation time 


qth scatter point 


Round-trip propagation time = tp/2-l-(tpq)-l-tq/2= (tp-l-tq)/2 -I- tpq 

Figure 12 Propagation Path of Indirect Scattering Response from Multiple 

Scatterers’ Interactions 
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Using equation 3.12, this second component can also possibly be modeled 
mathematically as a summation of the reflected echo signals from all of the multiple 
scatterers’ interactions: 

m m ~i 

^2(®)=E 

p=\ q=l,qjtp 

+ ...higher - order interaction 

where 

Apq is the amount of coupling between the pih and ^th scattering points, 

y/q is the phase (uniformly distributed over In) of the of ^th scatter, 

tq is the round-trip propagation time between the radar and the ^th scattering 
centers and 

tpq is the interaction time between the p\h and ^th scattering centers. 

Since the interaction time, tpq between the p\h and ^th scattering centers is 
generally unknown for complex targets, it is assumed to be equivalent to (tp+tq)/2 for the 
present (although this can be easily generalized). The radar backscattered signal due to 
multiple interactions is then scattered at a round-trip propagation time of (tp+tq)/2-l-(tp+tq)/2 
= tp+tq and the second major component of this signal can be rewritten as: 


¥+¥+(o\ 




'-+t 


mm { W 

«2(®)=Z «.)1 

p^l q=\,q^p 

+.... higher - order interaction 


(3.15) 


3. Mathematical Model for Combined Direct Specular and Indirect 
Multiple Scattering Interactions 

Collectively, the coherent radar backscattered signal of the target, H(a)) at radial 
frequency, a)=27rf, is a combination of both the direct specular and multiple scatter 
interaction terms of Hi(a>) and H 2 (ct>). 

Ignoring the effect of higher-order interactions and assuming that the signal 

amplitudes are frequency-independent, the radar backscattered signal, H(co) is given as: 
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(3.16) 


H{co) 




p=\ 




+ 


z 

p=lq=l,q^p 


-i[y/+ip+co{t^+tJ 


pq 


However in target range profiling, the point scatterer data are often observed in 
terms of relative range instead of the round-trip propagation times. As such, the radar 
backscattered signal can be rewritten using the radar range relationship of t=2r/c as: 


H(f) 


m - 

\A7iR—+y/ \ m 

m -i 




'Tj^pq^ 

c 

p=i 

p=\q-- 

^l,q*p 



(3.17) 


where 


Rp is the range between the radar and the p\h scatter point, 

Rq is the range between the radar and the ^th scatter point, 

/ is the frequency and 

c is the speed of light. 


D. EXAMPLE: RADAR BACKSCATTERED SIGNAL FROM A SIMPLE 3- 
POINT SCATTERERS TARGET 

To illustrate mathematical modeling of the backscattered radar signal, consider a 
3-point scatterers hypothetical target with scatter points at ranges of Ri, R2 and R3 from 
the radar as shown in Figure 13. 



Range, R 3 


Figure 13 Direct and Indirect Propagation Path of a 3-Point Scatterers Target 


24 


















Using equation 3.17, the radar backscatter signal for a given frequeney,/is given 


as: 


H(f) = 


3 -/ AttR —+y/ I ^ 3 


Z Z 


4fr{R^+R^)Uv'p+¥, 


p=l 


A^e 


f 


-i\ 4;rRi—+i//^ 


+ A^2^ 


+ A® 

f 

-i| A!i(Ri+R2)—+y/i+y/2 


p=\ q=\,q*p 
+1^2 


+ A^b 


■|42r(«| 


+ ^42^6 


f 

-('I 4n:(R^+R2)—+l|/^+l|/2 




f 

—zj ^— vy/^+y/-^ 




■ A^^b 


A 22 B 


f 

-i\ 471 R2—+y/2 
c 

f 

4a'(-Si +^3)—+((/] +(^3 
f 

47r(R2 +R2 )—+V'2 +¥2 
j^4»-(ff2 +^3 )^+!^2+!^3 


When 7?/ = 0 m, 7?2 = 1.5 m and Rs = 3.6 m, the corresponding radar backscatter 
signal can be rewritten as: 


H(f) =C, + Cf 


i\ 47l(\.5)^- 

-i 

4.>r(3.6)C 

-z|4.2-(5.1)^ 

L 

+ C-,B 

’ c _ 

+ C.B 


where 


Cl, Cl, C 3 and C 4 are constants. 


E. EXAMPLE: RADIAL RANGE PROFILE OF A SIMPLE 3-POINT 

SCATTERERS TARGET 

To better understand and visualize the contributions of individual scatterers and 
their multiple interactions in the radar backscattered signal, the radial range profile of the 
target is examined. The radial range profile of any target can be generated by evaluating 
the amplitudes and phases of the radar backscattered signals over a frequency band of 
evenly spaced frequency intervals and performing Fourier Transform to convert the 
signals into the time domain (Jouny, 1994). The radial range profile is a plot of the 
magnitude of the impulse response versus the projected down range along the line-of- 
sight between the radar system and the target, where scatterings from the individual 
specular scatterers and their multiple interactions will appear as peak responses. 


25 



The radial range profile of the simple 3-point seatterers target, computed using 
MATLAB’s Fast Fourier Transform of the associated radar backscattered signals over a 
frequency band from 5.0 GHz to 6.0 GHz, is illustrated in Figure 14. The related 
MATLAB program, “backscatter_signal.m” and source code used to generate this plot 
can be found in Appendix A. 



Figure 14 Radial Range Profile of 3-Point Seatterers Target 


From Figure 14, it can be observed that the first three peak impulse responses are 
consistent with the relative positions and geometry of the 3-points seatterers at ranges Ri 
= 0 m, ^2 = 1.5 m and R 3 = 3.6 m, while the fourth peak impulse response is a “ghosf’ 
artifact that does not map to any seatterers and represents the collective summation of the 
multiple scatters interaction terms between the 2"“^ and 3"^‘* scattering points at a range 
equivalent to R 2 + ^ 5 = 1.5 m + 3.6 m = 5.1 m. (Borden, 1998; Jouny, 1994). 
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IV. BI-SPECTRAL NCTR TECHNIQUE 


A. OVERVIEW OF BI-SPECTAL METHOD 

While the radial range profile NCTR teehnique proves to be a viable means of 
performing target reeognition, its applieation is predominantly based on the eorreet 
mapping and eorrelation of the peak impulse responses from the direet speeular seatterers 
to the aetual physieal layout of the target strueture and its eharaeteristie features. 

However, as diseussed and observed in the last seetion of Chapter III, some of 
these peak impulse responses in the radial range profile do not eorrespond to any target 
struetures and represent the multiple interaetions between major individual seatterers 
instead. Thus, a bi-speetral method, to differentiate the direet speeular seattering 
responses from the multiple seatterers interaetion responses, is proposed in Jouny (1994), 
Botha & Spoelstra (1992), Garber, Jouny & Moses (1995) and Jounty & Walton (1990) 
as another HRR imaging elass of NCTR teehnique to augment the radial range profile in 
performing radar target reeognition and identifieation. 

1. Definition of Bi-Spectrum 

In Higher-Orders Statistie (HOS) signal proeessing, the bi-speetrum of a 
eontinuous-time signal, h{\) is defined as the Fourier Transform of the third-order 
eumulant-generating funetion and is given in Jouny (1994), Botha & Spoelstra (1992), 
Garber, Jouny & Moses (1995) and Jounty & Walton (1990): 

00 00 

— 00 — CO 

where 

C 3 represents the triple eorrelation of the signal and is expressed as: 

C^{h,t2)= Ih"" {t)h{t + ty)h{t + t2)dt = Ef* {t)h{t + ti)h{t + t2)] 

— GO 

E{} denoted the expeetation and 

* superseript asterisk denoted a eomplex eonjugate. 
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2. Definition of Bi-Time 

However, in radar signal proeessing, the radar baekseattered signals are often 


proeessed and analyzed as diserete signals in the frequeney-domain (i.e. expressed as a 
function of frequency). The third-order cumulant or triple correlation of the frequency- 
domain signal, H{o)) is then defined by Jouny (1994), Botha & Spoelstra (1992), Garber, 
Jouny & Moses (1995) and Jounty & Walton (1990) as: 

00 

C 3 , 02 ) = j H*{a))H{a) + coi )H{a) + 02 )do) 

—00 (4-3) 

+ (yi)//(o + (»2)} 

The equivalent bi-spectrum, expressed as a function of time after applying the 
Fourier Transform, is known as the bi-time. The bi-time of a discrete time-domain bi- 
spectral signature is then expressed as: 

00 00 , . 

B{ti,t2)= X (4.4) 

0 )^ =-00 0 )^ =-00 

and can be represented in terms of impulse response as: 

B[ti,t2)^ E^[h)Kh)h*ih +G)} (4-5) 

where 

h{\) is the impulse response as a function of time 


3. Definition of Bi-Range 

The bi-time can be further expressed as a function of range using the radar range 
and round-trip propagation time relation of: 


r = 


ct 

~2 




t 



(4.6) 


where 

r is the range from the radar to the target, 

t is the round-trip propagation time needed for the radar signal to propagate 

to the target and back, and 

28 



c 


is the speed of light. 


Substituting t=lrlz and / =cd2n into the expression for bi-time, derives the 
eorresponding bi-range as (Jouny (1994), Botha & Spoelstra (1992), Garber, Jouny & 
Moses (1995) and Jounty & Walton (1990): 
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00 00 
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(4.7) 


and bi-range ean then be represented in terms of impulse response as: 

5(ri, r2) = E\i{ri )h{r2)h* {ri+r2)] (4.8) 

where 


h{x) is the impulse response as a funetion of range. 

From equation 4.7, it ean be dedueed that if the range profde for a given radar 
target eontains three responses, h{r\), hiri) and h{r^), two of whieh, h{r\) and h{r 2 ) are 
due to direet speeular seattering, and one is due to an interaetion between the direet 
seatterers, h{r^) , the bi-range will only deteet this term if the interaetion term appears at a 
distanee from the radar equal to the sum of the distanees of the direet terms, i.e. at rj, = ri 
+ ri aeeording to Jouny (1994). 


B. EXAMPLE: BI-RANGE OF A SIMPLE 3-POINT SCATTERERS TARGET 

Using the same simple 3-point seatterers target example as deseribed in Chapter 
III, the bi-range of the target ean be eomputed mathematieally using equation 4.7 as: 

B{ri , r2) = E\i{ri )h{r2 )h* {ri + r2 )} 

^ 1 / \ ( \] (4.9) 

~ ^ ^p^q^pq 1 “ ^p ’ ^2 “ ^q j"*" ~ ^q ’ ^2 “ ^p j. 

p^lq^l 
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Expanding the terms for the above expression: 


+ ^1^2^12 [^('l “ ^15 ^2 “ ^2 ) + ~ R2^r2 - R \)] 

+ A^sAs [^('l “ ^1 ’ ^^2 “ ^3 ) + ^('l “ ^3 ’ ^^2 “ ^1)] 

+ ^1^4^14[^('l ~R \^^2 “^4) + ^('1 “^4’''2 “^l)] 

+ ^2^1^21 [^('l “ ^2 ’ ^^2 “ ^1) + ^('l ~ R\^^2 ~ ^2 )] 

+ ^2^2^22 [^('l “ ^2 ’ ^^2 “ ^2 ) + ^('l “ ^2 ’ ^^2 “ ^2 )] 

+ ^2^3^423 [^('l ~R 2-'^2 “^3)+^('1 “^3’''2 “^ 2 )] 

+ ^^2^4^24 [^('l “ ^2 ’^^2 “ ^4 ) + ^('l ~ Ra^^ 2 ~ ^2 )] 

+ AAA\[^{^\ ~R3-'^2 “^i)+^('1 ~RI '>^2 “^ 3 )] 

+ ^ 3 ^ 2^432 [^('l ~R3-'^2 “^2)+^('1 “^2’'2 “^ 3 )] 

+ AAA3[^{’'1 ~R3^^2 “^ 3 ) + ^('1 ~R3^^2 “^ 3 )] 

+ ^ 3 ^ 4^34 [^(^1 -R3^^2 -Ra)+^[^\ ~Ra^^2 “^ 3 )] 

+ A^AiA^i \p{^\ ~ Ra^^2~ R\)^ ^('l “ ^1, ^2 “ ^4 )] 

+ ^4^2^42[^('l “^4’'2 “^2) + ^('l “^2’^2 “^ 4 )] 

+ ^4^3^43 [^(^1 -R 4 ^f ~2 “^3)+^('l ~R 3^^2 “^4)] 

+ A^A^A^4\d{ri -i?4,r2 -R^) + d{ri -R^,r2 -f?4)] 

Substituting Rj = 0 m, R 2 = 1.5 m, Rs = 3.6 m and R 4 =5.1 m, the bi-range ean 
be rewritten as: 

B{n.'■l) = C{[S(r^,r^)] + C2•[S(r^ -1.5,rj -l.S)] 

+ <^3'Kn -3.6,r 2 -3.6)] + C 4 '[^(ri -5.1,r2 -5.1)] 

+ C5'Kri,r2-1.5) + ^(ri-1.5,r2)] 

+ C^'[d{ri,r 2 -3.6) + ^(ri -3.6,r2)] 

+ C7'Kri,r2-5.l) + ^(ri-5.1,r2)] 

+ C^'[d[ri-\.5,r2 -3.6) + ^(ri -3.6,r2 -1.5)] 

+ C()'\d{ri -1.5,r2 -5.l)-l-(5’(rj -5.1,r2 -1.5)] 

+ Qo'Kn -3.6,r2 -5.l) + ^(ri -5.1,r2 -3.6)] 

where Ci’, C2’, C3’, C4’, C5’, Ce’, C7’, Cg’, C9’ and Cio’ are constants. 

Graphically, the bi-range of the 3-point scatterers target can be computed using 

MATLAB’s two-dimensional (2-D) discrete Fourier Transform of the third-order 

cumulant function of the radar backscatter signal. The results are then plotted as the 

three-dimensional (3-D) and 2-D (contour plot) bi-range profiles of the 3-point scatterers 
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target in Figure 15 and 16 respeetively. The related MATLAB programs, “bispeetral.m” 
and souree eode used to generate these plots ean be found in Appendix B. 



Figure 15 3-D Bi-Range Profile of 3-Point Seatterers Target 



Figure 16 2-D Bi-Range Profile of 3-Point Seatterers Target 


31 

















































































From Figures 15 and 16, it can be observed that the bi-range profiles are 
generated by responses at the range pairs (i?p and Rq) for p and ^ = 1,.... 4 with amplitude 
of ApAqApq. The locations of the responses are also noted to be consistent with the 
theoretically computed positions. 

C. CHARACTERISTIC BI-RANGE PROFILE 

Since a non-zero bi-range response is an outcome arising from the correlated 
interaction between the responses at range n and r 2 , it is clear that the peak responses in 
the bi-range profile of the radar backscattered signals represent strong additive multiple 
interactions between the scatterers or subcomponents of the target. 

Therefore, the characteristic bi-range profile of a target, which highlights the 
defining peak responses (within an order of magnitude) at coordinates (ri, r^) can 
potentially be employed to characterize the unique composition and layout of the target 
for recognition and identification. The characteristic bi-range profile, showcasing the 
peak responses of the 3-points-scatterer target example is illustrated in Figure 17. 



Figure 17 2-D Characteristic Bi-Range Profile of 3-Point Scatterers Target 
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It can be observed from Figure 17 that the peak responses in the bi-range profile 
of the 3-points seatterer target at {r\, ri) = (1.5m, 3.6m) eorresponded to the resulting 
multiple interactions between the 2"^^ scattering point at 1.5m and 3'^‘* seattering point at 
3.6m, as expeeted. 

The merits and limitations of bi-range profile NCTR technique aecording to Bao, 
Shi & Zhang (2001) and Bao, Pei & Xing (2001) are summarized in Table 4: 


Advantages 

Disadvantages 

a. Retains both the amplitude and phase 

information of a signal. Hence bi-spectral 

signatures are translation (or shift) invariant and 

robust to aspect dependence. 

a. Requires intensive computation of two- 

dimensional high range resolution data. The 

b. Good clutter and noise rejection 

performance due to ability to statistically suppress 

the clutter or noise characterized by zero-mean 

Gaussian probability density distribution. 

number of bi-spectral data is the square of the HRR 

profde. 

Tabled Summary of Bi-range Profi 

le’s Advantages and Disadvantages 
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V. BI-SPECTRAL NCTR EXPERIMENTATION 


A. “BACKHOE DATA DOME” EXPERIMENTATION 

To further examine the potential of using the eharaeteristie bi-range profile of a 
complex target for NCTR and identification, experimentation based on the Air Force 
Research Laboratory (AFRL), “Backhoe Data Dome, Version 1.0” by Lin & Naidu 
(2004) are performed and presented in the following sections. 

I. AFRL “Backhoe Data Dome Version I.O” 

The “Backhoe Data Dome, Version 1.0” consists of simulated wideband (7-13 
GHz), full polarization, complex backscatter data from a backhoe loader generated using 
a Computer-Aided-Design (CAD) model in free space as shown in Figure 18. The 
MATLAB program “plotfacet.m” available from the Air Force Research Laboratory 
Sensor Data Management System was used to plot all 3-D graphics of the backhoe 
loader. 



Figure 18 CAD Model of Backhoe Loader 


The backscatter data are then synthetically generated over a full 27i steradian 
viewing hemisphere with 14 samples per degree in both azimuth and elevation. There is 
one sample every 11.75 MHz over the 6 GHz bandwidth and full polarizations; i.e. 
horizontal linear polarization on transmit and receive (HH), vertical linear polarization on 
transmit and receive (VV) and horizontal linear polarization on transmit and vertical 
linear polarization on receive (HV). 
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A representation of the k-space data dome is as shown in Figure 19. 
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Figure 19 Backhoe Data Dome Representation in k-Space (Lin & Naidu, 2004) 

In this experiment, two sets of 2-D backhoe backscattered data which contained 
VV, HH and HV polarized k-space data from a 110 degree azimuth cut (from 350 
degrees to 100 degrees in azimuth) available through Lin & Naidu (2004) are used, one at 
0 degree elevation and the other at 30 degree elevation. 

With the method described in the preceding Chapters, characteristic bi-range 
profiles analysis on the 2-D k-space data is performed for three different test cases, with 
the backhoe loader setup in the following orientation and using the specified polarization; 

a. Test Case 1: Orientation of backhoe loader at 0 degree azimuth and 0 
degree elevation with W, TUT and HV polarizations. 

b. Test Case 2: Orientation of backhoe loader at 90 degree azimuth and 0 
degree elevation with W, HH and HV polarizations. 

c. Test Case 3: Orientation of bac kh oe loader at 60 degree azimuth and 30 
degree elevation with VV, HH and HV polarizations. 
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2. Test Case 1: Setup and Results (0azimuth=O°, (peievation=0°) 

In Test case 1, the backhoe loader is orientated at 0 degree azimuth and 0 degree 
elevation with reference to the radar direction as depicted in Figure 20; 



Figure 20 Orientation of CAD-Modeled Backhoe Loader at Aspect Viewing 
Angles of Azimuth, (0azimuth)=O° and Elevation, ((peievation)=0° 

The radial range profiles of the backhoe loader are computed using MATLAB’s 
Fast Fourier Transform of the associated radar backscattered signals, while the 3-D and 2 
-D bi-range profiles are computed using MATLAB’s two-dimensional (2-D) discrete 
Fourier Transform of the third-order cumulant function of the radar backscatter signals. 

Both radial range and bi-range profiles are processed over a frequency bandwidth 
of 5.9058 GFlz (from 7.0472 to 12.953 GFlz) with VV, TUT and FIV polarizations for the 
target’s aspect viewing angle of azimuth, (Oazimuth) = 0° and elevation, (cpeievation) = 0°. 
The related MATLAB program, “test_caseLm” and source code used to generate this 
plot can be found in Appendix C. 
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Magnitude of Impulse Response Magnitude of Impulse Response 


a. Radial Range Profile (0azimuth=O ° g?eievation=0 °) 

The radial range profiles of the backhoe loader with VV, HH and HV 
polarizations are illustrated in Figure 21, 22 and 23 respectively. 



Figure 21 Radial Range Profile of Backhoe Loader at Aspect Angles of 
0azimuth=O° and (peievation=0° with W Polarizafion 



Figure 22 Radial Range Profile of Backhoe Loader at Aspect Angles of 
0azimuth=O° and (Peievation=0° with HH Polarization 
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Figure 23 Radial Range Profile of Backhoe Loader at Aspect Angles of 
0azimuth=O° and (peievation=0° with HV Polarization 


b, 3-D Bi-Rcing€ Profile (0azimuth~O J (Pelevation—0 ^ 

The corresponding 3-D bi-range profiles of the backhoe loader with W, 
and HV polarizations are shown in Figure 23, 24 and 25 respectively. 
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Figure 24 3-D Bi-Range Profile of Backhoe Loader at Aspect Angles of 

0azimuth=O° and (Peievation=0° with VV Polarization 
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Figure 25 3-D Bi-Range Profile of Baekhoe Loader at Aspeet Angles of 

0azimuth=O° and (peievation=0° witli HH Polarization 
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Figure 26 3-D Bi-Range Profile of Baekhoe Loader at Aspeet Angles of 

0azimuth=O° and (peievation=0° with HV Polarizafion 
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c. 2-D Characteristic Bi-Range Profile ( 0azimuth=O ° q>eievation=0 °) 

The resulting 2-D characteristic bi-range profiles of the backhoe loader 
with VV, HH and HV polarizations are shown in Figure 25, 26 and 27 respectively. 



Figure 27 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth^O° and (peievation=0° with VV Polarization 




Figure 28 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth^O° and (peievation^O° with HH Polarization 
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Figure 29 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth=O° and (peievation=0° with HV Polarization 


3. Test Case 2: Setup and Results (0azimuth=9O°, (peievation=0°) 

In Test Case 2, the backhoe loader is orientated at 90 degree azimuth and 0 degree 
elevation with reference to the radar direction as depicted in Figure 30; 



Figure 30 Orientation of CAD-Modeled Backhoe Loader at Aspect Viewing 
Angles of Azimuth, (0azimuth)=9O° and Elevation, ((peievation)=0° 
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The radial range profiles of the baekhoe loader are eomputed using MATLAB’s 
Fast Fourier Transform of the associated radar backscattered signals, while the 3-D and 2 
-D bi-range profiles are computed using MATLAB’s two-dimensional (2-D) discrete 
Fourier Transform of the third-order cumulant function of the radar backscatter signals. 

Both radial range and bi-range profiles are processed over a frequency bandwidth 
of 5.9058 GHz (from 7.0472 to 12.953 GHz) with VV, HH and HV polarizations for the 
target’s aspect viewing angle of azimuth, (Oazimuth) = 90° and elevation, (cpeievation) = 0°. 
The related MATLAB program, “test_case2.m” and source code used to generate this 
plot can be found in Appendix D. 


a. Radial Range Profile ( 0azimuth=9O ° (Pelevation~ 0°) 

The radial range profiles of the baekhoe loader with VV, HH and HV 
polarizations are illustrated in Figure 31, 32 and 33 respectively. 



Figure 31 Radial Range Profile of Baekhoe Loader at Aspect Angles of 
0azimuth=9O° and (peievation=0° with VV Polarization 
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Figure 32 Radial Range Profile of Baekhoe Loader at Aspeet Angles of 
0azimuth^9O° and (peievation^O° witli HH Polarization 



Figure 33 Radial Range Profile of Baekhoe Loader at Aspect Angles of 
0azimuth=9O° and (peievation=0° with HV Polarization 
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Magnitude of Bi-range Response 


b. 3-D Bi-Range Profile ( 0azimuth=9O ° ^elevation~0 ^ 

The corresponding 3-D bi-range profiles of the backhoe loader with W, 
HH and HV polarizations are shown in Figure 34, 35 and 36 respectively. 
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Figure 34 3-D Bi-Range Profile of Backhoe Loader at Aspect Angles of 

0azimuth=9O° and (peievation=0° with W Polarization 
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Figure 35 3-D Bi-Range Profile of Backhoe Loader at Aspect Angles of 

0azimuth=9O° and (peievation=0° with HH Polarization 
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Figure 36 3-D Bi-Range Profile of Baekhoe Loader at Aspeet Angles of 

0azimuth=9O° and (peievation=0° witli HV Polarization 


c. 2-D Characteristic Bi-Range Profile (0azimuth=9O ° (peievation=0 °) 

The resulting 2-D eharaeteristie bi-range profiles of the baekhoe loader 
with VV, HH and HV polarizations are shown in Figure 37, 38 and 39 respeetively. 
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Figure 37 2-D Characteristie Bi-Range Profile of Baekhoe Loader at Aspeet 

Angles of 0azimuth=9O° and (peievation=0° with VV Polarization 
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Figure 38 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth=9O° and (peievation=0° with HFl Polarization 



Figure 39 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth=9O° and (peievation=0° with HV Polarization 
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4. Test Case 3: Setup and Results (0azimuth=6O°, (peievation=30°) 

In Test Case 3, the baekhoe loader is orientated at 60 degree azimuth and 30 
degree elevation with reference to the radar direction as depicted in Figure 40; 



Figure 40 Orientation of CAD-Modeled Baekhoe Loader at Aspect Viewing 
Angles of Azimuth, (0azimuth)=6O° and Elevation, ((peievation)=30° 

The radial range profiles of the baekhoe loader are computed using MATLAB’s 
Fast Fourier Transform of the associated radar backscattered signals, while the 3-D and 2 
-D bi-range profiles are computed using MATLAB’s two-dimensional (2-D) discrete 
Fourier Transform of the third-order cumulant function of the radar backscatter signals. 

Both radial range and bi-range profdes are processed over a frequency bandwidth 
of 5.9058 GHz (from 7.0472 to 12.953 GHz) with VV, HH and HV polarizations for the 
target’s aspect viewing angle of azimuth, (Oazimuth) = 60° and elevation, (cpeievation) = 30°. 
The related MATLAB program, “test_case3.m” and source code used to generate this 
plot can be found in Appendix E. 
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a. Radial Range Profile ( 6azimuth=60 ° ^elevation~30 °) 

The radial range profiles of the backhoe loader with VV, HH and HV 
polarizations are illustrated in Figure 41, 42 and 43 respectively. 



Figure 41 Radial Range Profile of Backhoe Loader at Aspect Angles of 
0azimuth=6O° and (peievation=30° with W Polarization 



Figure 42 Radial Range Profile of Backhoe Loader at Aspect Angles of 
0azimuth=6O° and (peievation=30° with FlFl Polarization 
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Figure 43 Radial Range Profile of Baekhoe Loader at Aspeet Angles of 
0azimuth=6O° and (peievation=30° witli FIV Polarization 


b. 3-D Bi-Range Profile ( 0azimuth=6O ° ^elevation— 30 °) 

The corresponding 3-D bi-range profiles of the baekhoe loader with W, 
and HV polarizations are shown in Figure 44, 45 and 46 respectively. 



Bi-range, (m) 


Bi-range, (m) 


Figure 44 3-D Bi-Range Profile of Baekhoe Loader at Aspect Angles of 

0azimuth=6O° and (peievation=30° with W Polarization 
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Figure 45 3-D Bi-Range Profile of Baekhoe Loader at Aspeet Angles of 

0azimuth=6O° and (peievation=30° with HH Polarization 



Bi-range, r, (m) 

Figure 46 3-D Bi-Range Profile of Baekhoe Loader at Aspeet Angles of 

0azimuth=6O° and (peievation=30° with FIV Polarization 
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c. 2-D Characteristic Bi-Range Profile ( 0azimuth=9O ° peievation=0 °) 
The resulting 2-D characteristic bi-range profiles of the backhoe loader 
with VV, HH and HV polarizations are shown in Figure 47, 48 and 49 respectively. 



Figure 47 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth^6O° and (Peievation=30° with W Polarization 



Figure 48 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth=6O° and (Peievation=30° with HH Polarization 
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Figure 49 2-D Characteristic Bi-Range Profile of Backhoe Loader at Aspect 

Angles of 0azimuth=6O° and (Peievation=30° with FIV Polarization 

B. ANALYSIS OF “BACKHOE DATA DOME” EXPERIMENTATION 

In studying and comparing the various radial range and bi-range profiles 
computed using the “backhoe data dome” experimentation for Test Cases 1, 2 and 3, 
observations and analysis relating to the uses of radial range profiles, polarization effects 
and bi-range profiles for target recognition can be rationalized and are discussed in the 
following sections. 

I. Radial Range Profiles 

The magnitude of the impulse response versus the range plots of the 
experimentation results provide a good representation of the one-dimensional radial 
profile of the backhoe loader from different aspect angles. It can be observed that 
locations of the peak responses appear to coincide and correspond to the relative 
positions of dominant individual scattering centers or subcomponents, such as the 
operator cab, engine compartment, and loader bucket, etc., that present a large RCS and 
consequently a strong backscattered signal to the radar system as illustrated in Figure 50. 
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Figure 50 Comparison of Radar Range Profile with Actual Backhoe Loader 

While it may seem that the radial range profile can possibly provide a measure of 
the length dimension of a target of interest along its line-of-sight from the radar system, 
this inference is usually erroneous and not equivalent to the actual physical target’s 
length due to two effects, as mentioned in Tait (2005). Instead of directly reflecting the 
radar signal back toward the radar system, multiple scatterers comprising the target and 
their interactions can result in indirect reflection and delay the round-trip propagation 
time, which then produces a radial range profile that is longer than the target. 
Conversely, the other effect results in a shorter radial range profile compared to the actual 
physical target’s length. This disparity is due to the shadowing or masking of one part 
the target by another part. When a certain part of a target is masked at a particular aspect 
angle, the obscuration reduces the electromagnetic energy that scattered off of it, and 
significantly reduces its contribution to the radial range profile. 
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For example, as also observed in Figure 50, the dipperstick of the backhoe loader 
is obscured by the engine compartment and operator cab, which then reduces its 
corresponding response in the radial range profile and results in an apparently shorter 
backhoe loader. 

2. Polarization Effects 

The polarization of a radar signal is defined as the direction of the electric E-field 
as the electromagnetic wave propagates. When a target of interest is illuminated by a 
linearly polarized signal, two echo signals are reflected back to the radar system; one 
with the same polarization as the transmitted signal (co-polarized signal) and the other 
with an orthogonal polarization (cross-polarized signal). The basis of using polarization 
effects for target recognition is to extract the differences in the radar backscattered 
signals from a target that is illuminated with different polarization. Although this method 
seems to be promising for discriminating simple geometrically shaped targets, such as 
spheres, comer reflectors and straight lengths of wire, etc, it is not as forthcoming for 
practical complex targets, such as aircraft and vehicles. 

For example, it is noted that the 3-D and 2-D bi-range profiles of the backhoe 
loader for the different test cases are not affected by the polarization of the radar signal as 
shown by the similarities between the different plots with VV and HH polarizations. 

This is probably due to the multiple scatterers or subcomponents interaction 
within the radar resolution cell that collectively resulted in a composite radar 
backscattered signal with a different and independent polarization as the transmitted 
signal. In addition, it is also observed from the 3-D and 2-D bi-range profiles with HV 
polarization that the peak responses from multiple scatterers’ interactions seems to be 
accentuated by the cross polarization effect. 

3. Characteristic Bi-Range Profile 

Most importantly, it is obvious from the experimentation results of the different 
test cases that there exists significant peak responses in the 3-D and 2-D bi-range profiles, 
implying strong additive multiple interactions between the scatterers or subcomponents 
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of the backhoe loader, that are characteristic of their bi-range shapes and locations. The 
positions of the multiple scattering interaction terms, relative to the direct specular 
scattering terms, correspond directly to the physical separation distances between the 
major individual scatterers of the backhoe loader. 

Therefore, these characteristic bi-range profiles can be drawn on, to highlight and 
map the unique multiple scattering interaction features and natures of the target which 
may then be used to compare and match against a reference database that contains bi¬ 
range signatures of different target types of potential interest and at all aspect angles of 
interest for target recognition and identification. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Among the detection systems under development in the DD(X) program is 

the Multi-Function Radar (MFR). The future MFR will provide Non- 

Cooperative Target Recognition (NCTR) and Kill Assessment. 

— Rear Admiral Philip Balisle, U.S. Navy (2002) 

A. SUMMING UP OF THESIS RESEARCH 

Radar-based NCTR will be a key capability and cardinal requirement in future 
radar systems equipped onboard all high-valued military combat platforms and assets, 
such as fighter aircraft and destroyers, etc. While there are many conventional 
techniques for the incorporation of NCTR into the next-generation radar systems, this 
thesis endeavors to expand on the knowledge and utility of the bi-spectral method for 
effective and robust target recognition and identification. 

The first objective of this thesis was to study and mathematically model the direct 
and indirect multiple interactions of scattering mechanisms of radar backscattered signal. 
The fundamentals of radar basic operation, electromagnetic wave propagation and 
scattering mechanisms of electromagnetic waves, when a radar signal intersects any 
target in its propagation path, were first reviewed. This then provided an insight into the 
relationships between the target’s radar signature and its physical size, geometry, shape, 
orientation and material composition, which were the basis in formulating the 
mathematical model of the radar backscattered signal. 

A mathematical model which represents the direct specular and indirect multiple 
scattering interactions of the backscattered radar signal was developed and subsequently 
verified using an example of a simple three-point scatterers target and the MATLAB 
program “backscatter signal.m” to plot the radial range profile. 

The second objective of this thesis was to develop a radar NCTR technique of 
using the bi-spectral signatures of backscattered radar signals for characterization of a 
target’s multiple interactions features. The definitions and mathematical derivations of 
the bi-spectrum, bi-time and bi-range of the radar backscattered signal were introduced. 


57 



together with the concept of using the bi-spectral signatures, also known as the 
characteristic bi-range profiles of backscattered radar signal as a complementary add-on 
to augment the radial range profile for target recognition and identification. 

The characteristic bi-range profile was then studied and verified through the same 
example of a simple 3-points scatterer target and the MATLAB program, “bispectral.m” 
was developed to compute and plot the bi-range profiles of the third-order cumulant or 
triple correlation of the frequency-domain radar backscattered signal. 

To further investigate the potential of using bi-spectra technique for radar target 
recognition, experimentation based on the AFRL, “Backhoe Data Dome” was conducted. 
The results and findings for three different test cases were analyzed and found to be 
promising. 

B. RECOMMENDATIONS FOR FUTURE WORK 

The novel idea of using the bi-spectral method for radar target recognition, based 
on the characterization of a target’s multiple scatterer interactions term, is a relatively 
recent development, considering that the concept was initially conceived and proposed as 
a basis of recognition by Jounty & Walton (1990) during the early nineties. As such, 
many aspects of this proposed bi-spectal NCTR technique can be further researched, 
developed and verified. The following sections recommend two areas for which future 
research and investigation can be considered. 

1. Experimentation Based on Actual Radar Signature Measurements of 
Real Complex Target 

For the current thesis, due to the restricted access to classified information, the 
experimentation to examine the potential of using the characteristic bi-range profile for 
target recognition and identification is based on the AFRL’s public release of the 
“Backhoe Data Dome, Version 1.0,” which consists of simulated wideband (7-13 GHz), 
full polarization, complex backscatter data from a backhoe loader generated through the 
use of a CAD model in free space (Lin & Naidu, 2004). While this set of data facilitated 
the verification of the concept of using a bi-spectral method for target recognition, the 
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experimentation laeks fidelity in the eonsideration of the radar bi-speetral signatures and 
seattering meehanisms from a real eomplex target in an environment surrounded with 
elutter and noise. 

Given proper seeurity elearanee, more eomprehensive experimentation, based on 
elassified data sets of aetual radar baekseattered signal measurements from real eomplex 
targets ean be examined and analyzed in future researeh to better understand the 
robustness of the bi-speetral method and its elutter and noise suppression performanee. 

2. Bi-spectral Waveform Design 

Another area of researeh interest that ean be further investigated and verified for 
future work is the topie of bi-speetral waveform design. The design of the radar 
waveform is erueial in the measurement of the radar baekseattered signals required to 
generate the bi-speetral profile. In order to generate bi-range profiles, whieh are 
eompatible with the referenee signature database of different target types of potential 
interest for target reeognition, the radar waveform needs to provide high range resolution 
measurements and extraet diseriminating multiple seatterers interaetion features of the 
target. 

While the experimentation results from the eurrent thesis seem to suggest that the 
use of a eross polarized signal ean possibly aeeentuate the peak bi-range profde responses 
from multiple seatterers’ interaetions, more in-depth investigation and studies have to be 
performed in order to understand and validate the polarization eharaeteristies of the 
multiple seatterers’ interaetions terms. Apart from the polarization of the radar signal, 
other radar waveform design parameters like the transmitted frequeney, pulse width, 
pulse repetition frequeney, frequeney bandwidth and modulation, ete., ean also be 
examined in future researeh work in order to establish an optimal bi-speetral waveform. 
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APPENDIX 


MATLAB SOURCE CODES 


A. SOURCE CODE FOR RADAR BACKSCATTERED SIGNAL 

% - 

% backscattersignal.m 

% 

% This program computes the radar hackscattered signal and 
% generates the radial range profde of n points scatterers 

% 

% Written hy Yeo, Jiunn Wah 

% - 

clear all 


n = 3; 

% 

Number of point scatter 

f= 5:0.01:6; 

% 

Frequency band in GHz 

A=[l 1 1]; 

% 

Amplitude of the scattered signal 

R=[0 1.5 3.6]; 

% 

Range between the Radar and point scatter 

s = 0; 

% 

Initialization of counter 

c = 3e8; 

% 

Speed of light 


bandwidth = f(length(f))-f(l); 

% Direct specular scattering term 
for countl=l:n 

s=s+A(countl).*exp(i.*2.*pi.*2*R(countl).*le9.*f./c);%-rand*2*pi); 

end 


% Indirect multiple scattering interaction term 
for count2=l :n 
for count3=l:n 
if count3~=count2 

s=s+A(count2)*A(count3).*exp(i.*2.*pi.*2*(R(count2)+R(count3)).*le9.*f./c); 

end 

end 

end 


% Plotting of Radial Range Profde 

z=fft(s) ./length(s); 

figure 

plot(linspace(-length(z)/2,length(z)/2,length(z))*c/2/bandwidth/le9, fftshift(abs(z))); 
xlabel('Range (m)') 

ylabel('Magnitude Impulse Response') 
grid on 
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B. 


SOURCE CODE FOR GENERATING THE THIRD-ORDER CUMULANT 
AND BI-SPECTRAL PLOTS 


% - 

% bispectral.m 

% 

% This program generates the 3-D and 2-D bi-range profile of 

% a radar backscattered signal 

% 

% Written by Yeo, Jiunn Wah 

% - 

clear all 


% Call program to compute the radar backscattered signal and 
% generates the radial range profile of n points scatterer 
backscattersignal; 

% Function to generate the third order cumulant of the radar backscattered signal 
[ex] = cum3(s,(length(s)-l)); 

% Perform Fourier Transform of the third order cumulant to generate the bispectrum 
bisp = fftshift(fft2(ifftshift(cx))); 

% Setting the plot limits 

X = linspace(-length(bisp)/4,length(bisp)/4,length(bisp))*c/2/bandwidth/le9; 
y = hnspace(-length(bisp)/4,length(bisp)/4,length(bisp))*c/2/bandwidth/le9; 
ref = find(x==0); 

% Plotting of the 3-D bi-range profile 
figure 

mesh(x(ref length(x)),y(reflength(y)),abs(bisp(ref length(bisp),(ref length(bisp))))); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 
zlabel('Magnitude of Bi-range Response') 

% Plotting of the 2-D Characteristic bi-range profile 
figure 

contourf(x(ref length(x)),y(ref length(y)),abs(bisp(ref length(bisp),(ref length(bisp)))), [2e4 1.2e5 1.4e5 

1.6e5 1.8e5 2e5]); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

grid on 

colorbar 


% - 

% cum3.m 

% 

% This function generates the third order cumulant of the radar backscattered signal 

% 

% Written by Yeo, Jiunn Wah 

% 

% ex : the cross cumulant matrix C(m,n) of the radar backscattered signal 

% signal : the radar backscattered signal to be processed 

% 

% signal length : the maximum argument of the cross cumulant 

% - 
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function [cx] = cum3 (signal, signal length) 

% Assignment and pre-conditioning of signal 
signal_ 1 =conj (signal(:)); 
signal_2=signal(:); 
signal_3=signal(:); 

N = lengtli(signal_l); 

% Triple correlation of signal 

XI = repmat(signal_l,l,2*signal_length-l-l); 

C2=[zeros(signal_lengtli, 1); signal_2( 1: signallength-l-1)]; 
R2=[signal_2(signal_length-l-l:N).' zeros(l,signal_length)]; 
X2=liankel(C2,R2); 

C3=[zeros(signal_lengtli, 1); signal_3 (1: signallength-l-1)]; 
R3=[signal_3(signal_length-l-l:N).' zeros(l,signal_length)]; 
X3=hankel(C3,R3).'; 

cx=X2*(Xl.*X3)/N; 


C. SOURCE CODE FOR TEST CASE 1 

% - 

% test_casel.m 

% 

% This program computes the third order cumulant and bispectrum of the “backhoe 
% data dome” radar backscattered signals for a target orientation of azimuth angle, 

% (0azimuth) = 0° and elevation anlge, ((Peievation) = 0°. 

% 

% Both radial range and bi-range profiles are processed and plotted over a frequency 
% bandwidth of 5.9058 GHz (from 7.0472 to 12.953 GHz) with W, HH and HV 
% polarizations. 

% 

% Written by Yeo, Jiunn Wah 

% - 

clear all; 

% Load “Backhoe Data Dome” data 

load('C:\DOCUME~l\J&C\MYDOCU~l\NPS\Thesis\BACKHO~l\2D_CHA~l\2D_K_S~l\backhoe_elOO 
0_az3 50to 100 .mat') 


c = 3e8; 

% 

Speed of light 


bandwidth = 5.9058; 

% 

Stimulated frequency bandwidth 


sW 0 = W(:,141); 
zW_0 = fft(sVV_0); 

% 

W polarized radar signal at Oazimuth = 0° and 

^elevation ~ 0 

sHH 0 = HH(:,141); 
zHHO = fft(sHH_0); 

% 

HH polarized radar signal at Oazimuth = 0° and 

^elevation “ 0 

sHV_0 = HV(:,141); 

% 

HV polarized radar signal at Oazimuth = 0° and 

^elevation “ 0 
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zHV_0 = fft(sHV_0); 


% Plotting of Radial Range Profile with W polarization 
figure 

plot((-length(zVV_0)/2: length(zW_0)/2-1 )*c/2/bandwidth/l e9, (abs(fftshift(zW_0)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with VV Polarization') 
grid on 


% Plotting of Radial Range Profile with HH polarization 
figure 

plot((-length(zHH_0)/2:length(zHH_0)/2-l)*e/2/bandwidth/le9, fftshift(abs((zHH_0)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with HH Polarization') 
grid on 


% Plotting of Radial Range Profile with HV polarization 

figure 

plot((-length(zHV_0)/2:length(zHV_0)/2-l)*c/2/bandwidth/le9, (abs(fftshift(zHV_0)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with HV Polarization') 
grid on 


% Computation of third order cumulant and bispeetrum 
[cxW_0] = eum3(sVV_0,(length(sW_0)-l)); 
bispWO = fftshift(fft2(ifftshift(cxW_0))); 

[cxHHO] = eum3(sHH_0,(length(sHH_0)-l)); 
bispHHO = fftshift(fft2(ifftshift(cxHH_0))); 

[cxHV_0] = eum3(sHV_0,(length(sHV_0)-l)); 
bispHVO = fftshift(fft2(ifftshift(cxHV_0))); 

% Plotting of the 3-D bi-range profile with W polarization 
figure 

xW_0=linspace(-length(bispW_0)/4,length(bispVV_0)/4,length(bispW_0))*c/2/bandwidth/le9; 
yW_0 = bnspace(-lengtIi(bispW_0)/4,lengtIi(bispVV_0)/4,lengtIi(bispW_0))*c/2/bandwidth/le9; 
refW_0 = fmd(xW_0==0); 
refendW_0 = find(xW_0>=4.5 & xW_0<=4.52); 

mesh(x W_0(refW_0 :refendW_0),yW_0(refW_0 :refendW_0),abs(bispW_0(refW_0 :refendW_0, 

(refW_0:refendW_0)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with W Polarization') 

% Plotting of the 3-D bi-range profile with HH polarization 
figure 

xHHO = hnspaee(-length(bispHH_0)/4,length(bispHH_0)/4,length(bispHH_0))*c/2/bandwidth/le9; 
yHHO = hnspaee(-length(bispHH_0)/4,length(bispHH_0)/4,length(bispHH_0))*c/2/bandwidth/le9; 
refHHO = fmd(xHH_0==0); 
refendHH O = find(xHH_0>=4.5 & xHH_0<=4.52); 
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mesh(xHH_0(refHH_0 :refendHH_0),yHH_0(refHH_0 :refendHH_0),abs(bispHH_0(refHH_0 :refendHH_0, 

(refHH_0:refendHH_0)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

tide('3-D Bi-range Profde for \theta_a_z = 0 deg and \theta_e_l = 0 deg with HH Polarization') 

% Plotting of the 3-D bi-range profde with HV polarization 
figure 

xHVO = linspace(-length(bispHV_0)/4,length(bispHV_0)/4,length(bispHV_0))*c/2/bandwidth/le9; 
yHVO = linspace(-length(bispHV_0)/4,length(bispHV_0)/4,length(bispHV_0))*c/2/bandwidth/le9; 
refHV_0 = find(xHV_0>=-2.51 & xHV_0<=-2.5); 
refendHV_0 = find(xHV_0>=2.5 & xHV_0<=2.51); 

mesh(xHV_0(refHV_0:refendHV_0),yHV_0(refHV_0:refendHV_0),abs(bispHV_0(refHV_0:refendHV_0, 

(refHV_0:refendHV_0)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with HV Polarization') 

% Plotting of the 2-D characteristic bi-range profile with VV polarization 
figure 

contour(xVV_0(refW_0:refendW_0),yW_0(refW_0:refendW_0),abs(bispW_0(refW_0:refendW 
_0,(refW_0:refendW_0))), [2e7 4e7 6e7 8e7 le8 1.2e8 1.4e8 1.6e8 1.8e8 2e8]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Characteristic Bi-range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with W 

Polarization') 

grid on 

colorbar 


% Plotting of the 2-D characteristic bi-range profile with HH polarization 
figure 

contour(xHH_0(refHH_0:refendHH_0),yHH_0(refHH_0:refendHH_0),abs(bispHH_0(refHH_0:refendHH 
_0,(refHH_0:refendHH_0))), [2e7 4e7 6e7 8e7 le8 1.2e8 1.4e8 1.6e8 1.8e8 2e8]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Characteristic Bi-range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with HH 

Polarization') 

grid on 

colorbar 


% Plotting of the 2-D characteristic bi-range profile with HV polarization 
figure 

contour(xHV_0(refHV_0:refendHV_0),yHV_0(refHV_0:refendHV_0),abs(bispHV_0(refHV_0:refendHV 
_0,(refHV_0:refendHV_0))), [3 5 10 15 20 25 30]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Characteristic Bi-range Profile for \theta_a_z = 0 deg and \theta_e_l = 0 deg with HV 

Polarization') 

grid on 

colorbar 
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D. SOURCE CODE FOR TEST CASE 2 


% - 

% test_case2.m 

% 

% This program computes the third order cumulant and bispectrum of the “backhoe 
% data dome” radar backscattered signals for a target orientation of azimuth angle, 

% (dazimuth) = 90° and elevation anlge, ((peievation) = 0°. 

% 

% Both radial range and bi-range profdes are processed and plotted over a frequency 
% bandwidth of 5.9058 GHz (from 7.0472 to 12.953 GHz) with W, HH and HV 
% polarizations. 

% 

% Written by Yeo, Jiunn Wah 

% - 

clear all; 


% Load “Backhoe Data Dome” data 

load('C:\DOCUME~l\J&C\MYDOCU~l\NPS\Thesis\BACKHO~l\2D_CHA~l\2D_K_S~l\backhoe_elOO 
0_az3 50to 100 .mat') 


c = 3e8; 

% 

Speed of light 


bandwidth = 5.9058; 

% 

Stimulated frequency bandwidth 


sW 90 = VV(:,1401); 
zVV_90 = fft(sVV_90); 

% 

W polarized radar signal at Oazimuth = 90° and 

^elevation “ 0 

sHH 90 = HH(:,1401); 
zHH_90 = fft(sHH_90); 

% 

HH polarized radar signal at Oazimuth = 90° and 

^elevation ~ ^ 

sHV 90 = HV(:,1401); 
zHV_90 = fft(sHV_90); 

% 

HV polarized radar signal at Oazimuth = 90° and 

^elevation ^ 0 


% Plotting of Radial Range Profde with W polarization 
figure 

plot((-length(zW_90)/2 :length(zW_90)/2-1 )*c/2/bandwidth/1 e9, (abs(fftshift(zW_90)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with W Polarization') 
grid on 


% Plotting of Radial Range Profile with HH polarization 
figure 

plot((-length(zHH_90)/2:length(zHH_90)/2-l)*c/2/bandwidth/le9, fftshift(abs((zHH_90)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with HH Polarization') 
grid on 


% Plotting of Radial Range Profile with HV polarization 
figure 

plot((-length(zHV_90)/2 :length(zHV_90)/2-1 )*c/2/bandwidth/1 e9, (abs(fftshift(zHV_90)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with HV Polarization') 
grid on 
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% Computation of third order cumulant and bispectrum 
[cxW_90] = cum3(sW_90,(length(sW_90)-l)); 
bispW_90 = fftshift(fft2(ifftshift(cxW_90))); 

[cxHH_90] = cum3(sHH_90,(length(sHH_90)-l)); 
bispHH_90 = fftshift(fft2(ifftshift(cxHH_90))); 

[cxHV_90] = cum3(sHV_90,(length(sHV_90)-l)); 
bispHV_90 = fftshift(fft2(ifftshift(cxHV_90))); 

% Plotting of the 3-D bi-range profde with W polarization 
figure 

xW_90 = linspace(-length(bispW_90)/4,length(bispW_90)/4,length(bispW_90))*c/2/bandwidth/l e9; 
yW_90 = linspace(-length(bispW_90)/4,length(bispW_90)/4,length(bispW_90))*c/2/bandwidth/l e9; 
refW_90 = fmd(xW_90==0); 
refendW_90 = find(xW_90>=2 & xW_90<=2.0085); 

mesh(xW_90(refW_90 :refendW_90),yW_90(refvV_90 :refendW_90),abs(bispW_90(refW_90 :ref 

endW_90,(refW_90:refendVV_90)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with W Polarization') 

% Plotting of the 3-D bi-range profile with HH polarization 
figure 

xHH_90 = linspace(-length(bispHH_90)/4,length(bispHH_90)/4,length(bispHH_90))*c/2/bandwidth/le9; 
yHH_90 = linspace(-length(bispHH_90)/4,length(bispHH_90)/4,length(bispHH_90))*c/2/bandwidth/le9; 
refHH_90 = find(xHH_90==0); 
refendHH_90 = find(xHH_90>=2 & xHH_90<=2.0085); 

mesh(xHH_90(refHH_90:refendHH_90),yHH_90(refHH_90:refendHH_90),abs(bispHH_90(refHH_90:ref 

endHH_90,(refHH_90:refendHH_90)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with HH Polarization') 

% Plotting of the 3-D bi-range profile with HV polarization 
figure 

xHV_90 = linspace(-length(bispHV_90)/4,length(bispHV_90)/4,length(bispHV_90))*c/2/bandwidth/le9; 
yHV_90 = linspace(-length(bispHV_90)/4,length(bispHV_90)/4,length(bispHV_90))*c/2/bandwidth/le9; 
refHV_90 = find(xHV_90>=-2.01 & xHV_90<=-2.0); 
refendHV_90 = find(xHV_90>=2.0 & xHV_90<=2.01); 

mesh(xHV_90(refHV_90 :refendHV_90),yHV_90(refHV_90 :refendHV_90),abs(bispHV_90(refHV_90 :ref 

endHV_90,(refHV_90:refendHV_90)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with HV Polarization') 

% Plotting of the 2-D characteristic bi-range profile with VV polarization 
figure 

contour(xVV_90(refW_90:refendW_90),yW_90(refW_90:refendW_90),abs(bispW_90(refW_90:r 
efendW_90,(refW_90:refendW_90))), [2.5e9 5e9 7.5e9 lelO 1.5el0 2el0 2.5el0]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 
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title('2-D Characteristic Bi-range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with VV 

Polarization') 

grid on 

colorbar 


% Plotting of the 2-D characteristic bi-range profile with HH polarization 
figure 

contour(xHH_90(refHH_90:refendHH_90),yHH_90(refHH_90:refendHH_90),abs(bispHH_90(refHH_90:r 
efendHH_90,(refHH_90:refendHH_90))), [2.5e9 5e9 7.5e9 lelO 1.5el0 2el0 2.5el0]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Characteristic Bi-range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with HH 

Polarization') 

grid on 

colorbar 


% Plotting of the 2-D characteristic bi-range profile with HV polarization 
figure 

contour(xHV_90(refHV_90 :refendHV_90),yHV_90(refHV_90 :refendHV_90),abs(bispHV_90(refHV_90 :r 
efendHV_90,(refHV_90:refendHV_90))), [18 20 40 60 80 100 120 140 160 180]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Characteristic Bi-range Profile for \theta_a_z = 90 deg and \theta_e_l = 0 deg with HV 

Polarization') 

grid on 

colorbar 

E. SOURCE CODE FOR TEST CASE 3 

% - 

% test_case3.m 

% 

% This program computes the third order cumulant and bispectrum of the “backhoe 

% data dome” radar backscattered signals for a target orientation of azimuth angle, 

% (fiazimuth) = 60° and elevation anlge, ((peievation) = 30°. 

% 

% Both radial range and bi-range profiles are processed and plotted over a frequency 
% bandwidth of 5.9058 GHz (from 7.0472 to 12.953 GHz) with W, HH and HV 
% polarizations. 

% 

% Written by Yeo, Jiunn Wah 

% - 

clear all; 


% Load “Backhoe Data Dome” data 

load('C:\DOCUME~l\J&C\MYDOCU~l\NPS\Thesis\BACKHO~l\2D_CHA~l\2D_K_S~l\backhoe_el03 
0_az350to 100 .mat') 


c = 3e8; 

% 

Speed of light 

bandwidth = 5.9058; 

% 

Stimulated frequency bandwidth 

sW 60 = VV(:,981); 
zVV_60 = fft(sVV_60); 

% 

W polarized radar signal at Oazimuth = 60° and (peievation = 30' 

sHH 60 = HH(:,981); 
zHH_60 = fft(sHH_60); 

% 

HH polarized radar signal at Oazimuth = 60° and (peievation = 30' 
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sHV_60 = HV(:,981); % 

zHV_60 = fft(sHV_60); 


HV polarized radar signal at Bazimuth = 60° and (peievation = 30' 


% Plotting of Radial Range Profde with W polarization 
figure 

plot((-length(zW_60)/2 :length(zW_60)/2-1 )*c/2/bandwidth/1 e9, (abs(fftshift(zW_60)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 60 deg and \theta_e_l = 30 deg with W Polarization') 
grid on 


% Plotting of Radial Range Profile with HH polarization 
figure 

plot((-length(zHH_60)/2:length(zHH_60)/2-l)*c/2/bandwidth/le9, fftshift(abs((zHH_60)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 60 deg and \theta_e_l = 30 deg with HH Polarization') 
grid on 


% Plotting of Radial Range Profile with HV polarization 
figure 

plot((-length(zHV_60)/2:length(zHV_60)/2-l)*e/2/bandwidth/le9, (abs(fftshift(zHV_60)))) 
xlabel('Range (m)') 

ylabel('Magnitude of Impulse Response') 

title('Radial Range Profile for \theta_a_z = 60 deg and \theta_e_l = 30 deg with HV Polarization') 
grid on 


% Computation of third order cumulant and bispectrum 
[cxVV_60] = cum3(sVV_60,(length(sVV_60)-l)); 
bispWhO = fftshift(fft2(ifftshift(cxVV_60))); 

[cxHH_60] = cum3(sHH_60,(length(sHH_60)-l)); 
bispHH_60 = fftshift(fft2(ifftshift(cxHH_60))); 

[cxHV_60] = cum3(sHV_60,(length(sHV_60)-l)); 
bispHV_60 = fftshift(fft2(ifftshift(cxHV_60))); 

% Plotting of the 3-D bi-range profile with W polarization 
figure 

xVV_60 = linspace(-length(bispVV_60)/4,length(bispVV_60)/4,length(bispVV_60))*c/2/bandwidth/l e9; 
yVV_60 = linspace(-length(bispVV_60)/4,length(bispVV_60)/4,length(bispVV_60))*c/2/bandwidth/l e9; 
refVV_60 = find(xVV_60==0); 
refendVV_60 = find(xVV_60>=2 & xVV_60<=2.01); 

mesh(xVV_60(refW_60:refendVV_60),yVV_60(refW_60:refendVV_60),abs(bispVV_60(refW_60:ref 

endVV_60,(refVV_60:refendVV_60)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \theta_a_z = 60 deg and \theta_e_l = 30 deg with VV Polarization') 

% Plotting of the 3-D bi-range profile with HH polarization 
figure 

xHH_60 = linspace(-length(bispHH_60)/4,length(bispHH_60)/4,length(bispHH_60))*c/2/bandwidth/le9; 
yHH_60 = linspace(-length(bispHH_60)/4,length(bispHH_60)/4,length(bispHH_60))*c/2/bandwidth/le9; 
refHH_60 = find(xHH_60==0); 
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refendHH_60 = fmd(xHH_60>=2 & xHH_60<=2.01); 

mesh(xHH_60(refHH_60:refendHH_60),yHH_60(refHH_60:refendHH_60),abs(bispHH_60(refHH_60:ref 

endHH_60,(refHH_60:refendHH_60)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

tide('3-D Bi-range Profde for \tbeta_a_z = 60 deg and \tbeta_e_l = 30 deg with HH Polarization') 

% Plotting of tbe 3-D bi-range profde with HV polarization 
figure 

xHV_60 = linspace(-lengtb(bispHV_60)/4,lengtb(bispHV_60)/4,lengtb(bispHV_60))*c/2/bandwidtb/le9; 
yHV_60 = linspace(-lengtb(bispHV_60)/4,lengtb(bispHV_60)/4,lengtb(bispHV_60))*c/2/bandwidtb/le9; 
refHV_60 = fmd(xHV_60>=-2.01 & xHV_60<=-2.0); 
refendHV_60 = find(xHV_60>=2.0 & xHV_60<=2.01); 

mesb(xHV_60(refHV_60:refendHV_60),yHV_60(refHV_60:refendHV_60),abs(bispHV_60(refHV_60:ref 

endHV_60,(refHV_60:refendHV_60)))); 

xlabel('Bi-range, r_l (m)') 

ylabel('Bi-range, r_2 (m)') 

zlabel('Magnitude of Bi-range Response') 

title('3-D Bi-range Profile for \tbeta_a_z = 60 deg and \tbeta_e_l = 30 deg with HV Polarization') 

% Plotting of tbe 2-D cbaracteristic bi-range profile with VV polarization 
figure 

contour(xVV_60(refW_60 :refendW_60),yW_60(refW_60 :refendW_60),abs(bispW_60(refW_60 :r 
efendW_60,(refW_60:refendW_60))), [1.8e5 4e5 6e5 8e5 le6 1.2e6 1.4e6 1.6e6 1.8e6]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Cbaracteristic Bi-range Profile for \tbeta_a_z = 60 deg and \tbeta_e_l = 30 deg with W 

Polarization') 

grid on 

colorbar 


% Plotting of tbe 2-D cbaracteristic bi-range profile with HH polarization 
figure 

contour(xHH_60(refHH_60:refendHH_60),yHH_60(refHH_60:refendHH_60),abs(bispHH_60(refHH_60:r 
efendHH_60,(refHH_60:refendHH_60))), [2.5e5 4e5 6e5 8e5 le6 1.2e6 1.4e6 1.6e6 1.8e6 2.0e6 2.5e6]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Cbaracteristic Bi-range Profile for \tbeta_a_z = 60 deg and \tbeta_e_l = 30 deg with HH 

Polarization') 

grid on 

colorbar 


% Plotting of tbe 2-D cbaracteristic bi-range profile with HV polarization 
figure 

contour(xHV_60(refHV_60 :refendHV_60),yHV_60(refHV_60 :refendHV_60),abs(bispHV_60(refHV_60 :r 
efendHV_60,(refHV_60:refendHV_60))), [3.5e6 4e6 6e6 8e6 le7 2e7 2.5e7 3e7 3.5e7]); 
xlabel('Bi-range, r_l (m)') 
ylabel('Bi-range, r_2 (m)') 

title('2-D Cbaracteristic Bi-range Profile for \tbeta_a_z = 60 deg and \tbeta_e_l = 30 deg with HV 

Polarization') 

grid on 

colorbar 
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